Half quantum vortex in superfluid ^He-A phase in parallel plate geometry 
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The half quantum vortex(HQV) in condensate has been studied, since it was predicted by Salomaa 
and Volovik in superfluid "^He-A phase. However, an experimental evidence for its existence has not 
been reported so far. Motivated by a recent experimental report by Yamashita et al^, we study 
the HQVs in superfluid '^He confined between two parallel plates with a gap D ~ 10 fj,m in the 
presence of a magnetic field H ~ 26 mT perpendicular to the parallel plates. We find that the 
bound HQVs are more stable than the singular vortices and free pairs of HQVs, when the rotation 
perpendicular to the parallel plates is below the critical speed, f2c ~ 2 rad/s. The bound pair of 
HQVs accompanies the tilting of d-vector out of the plane, which leads to an additional absorption 
in NMR spectra. Our study appears to describe the temperature and rotation dependence of the 
observed satellite NMR signal, which supports the existence of the HQVs in ''He. 

PACS numbers: 67.57.-z,67.57.Fg,67.57.Lm 



I. INTRODUCTION 

The superfluid "^He-A is a remarkable condensate with 
several broken symmetries and a variety of topologi- 
cal defects. It was proposed that many concepts and 
notions generated in the high energy particle physics 
can be tested in the earth bound ultra-low temperature 
laboratoriesi^i^. Indeed a variety of vortices, which 
include the continuous vortices (the vortices with soft- 
core) have been observed in the rotating superfluid '^He- 
All of the continuous vortices are associated with 
Z-texture are stable in the open systems where angular 
momentum of Cooper pair, ^-vector is not constrained. 

On the other hand, in the parallel plate geometry, 
where superfluid ^He is confined between two parallel 
plates with a gap D ^ ^u, where 10 ^m) is 

the dipolc coherence length, the /-vector is forced to be 
perpendicular to the plate.— In such a geometry, the 
continuous vortices cannot exist, and one may conclude 
that the vortices with the circulation = ±1 called 
the phase vortices or the singular vortices are only possi- 
ble vorticesi^ However, Salomaa and Voloviki^ predicted 
that the half quantum vortices(HQVs) with N — zt^ 
should exist in the parallel plate geometry. In spite of 
intensive experimental search, the HQV has not been 
found in the rotating superfluid '^He-A in the conflned 
geometry.— 

It is worthwhile to note that the HQVs have played 
an important role outside the superfluid '^Hc-A. For ex- 
ample, the HQVs in the tricrystal geometry in the high 
temperature cuprate superconductors played the crucial 
role in identifying d-wave symmetry of the high temper- 
ature superconducting order parameter.— The HQVs in 
triplet superconductors have been also suggeste d^^'^'^i^^ 
to interpret the extremely strong vortex pinning observed 
in triplet superconductors such as UPts, Ui-ajTh^jBeis, 
and Sr2Ru04'^ 

The present work is motivated by the recent experi- 
ment by Yamashita et alii We consider the superfluid 



■^He-A confined between the two parallel plate with the 
gap D ^ S.Dj where the rotation axis and the mag- 
netic field are perpendicular to the parallel plates. Us- 
ing the standard formula for free energies given in 
we found the foUowings. (a) The singular vortices with 
= ±1 are unstable and split into pairs of HQV 
with N = zt^, unless the rotation speed is very high, 
57 > 0„ ~ 5 X 10^ rad/s. (b) For small rotation speed 
(fi < ^--^ 2 rad/s); the bound pairs of HQVs arc more 
stable than the free HQVs. (c) Only the bound pairs of 
HQVs are detectable by NMR as satellite signals. Both 
the singular vortices and free HQVs lead to NMR sig- 
nals indistinguishable from backgrounds, since they do 
not disturb the configuration of / _L d. (d) Our results 
of the satellite frequencies and the intensities are consis- 
tent with the observation of satellite peak in the parallel 
plate geometry.— In particular, while flc weakly depends 
on the temperature and the magnetic field strength, both 
the satellite frequency and the satellite intensity are al- 
most independent of the field strength and temperature. 



II. FREE ENERGIES OF ONE SINGULAR 
VORTEX AND A PAIR OF HQV 



The superconducting order parameter of superfluid 
■^Hc-A is characterized by 2 unit vectors called 1-vector, 
I and d-vector, d. We can write the gap function 



where 



A,b(k) = D(k) • iaia2)ab 



D{k) = M{ki + ik2) 



(1) 



(2) 



Here (/.j = 1,2,3) are Pauli matrices and a, b are 
spin I and [. ki (i = 1,2) arc the projection of the unit 
wave vector k in the plane perpendicular to I. Therefore 
I corresponds to the quantization axis of the angular mo- 
mentum of the Cooper pair. This is also called the chiral 
vector. 
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When the spatial variation of d and I are much slower 
than ^ (~ 0.01/j,to), coherence length, the spatial con- 
figuration of I and d arc described by the texture free 
energy.—. In the following, we shall compare the texture 
free energy of a singular vortex(or phase vortex), and a 
free pair of HQV with an optimized distance R, where R 
is a distance between two HQVs. 

The texture free energy is given by 



F = \xnC^ I dxdy j i^(V$)2 + ^ \dd, 

■(3) 

where we assumed the parallel plates extend in the x — y 
plane. Here xn and C are spin susceptibility and spin 
wave velocity. Then K is given by the ratio between the 
superfluid density, ps and spin supcrfluid density, psp. 



K = 



Ps 
Psp 



1 



(4) 



where Fi and are Landau's Fermi liquid coefficients. 
We also assumed I \\ z everywhere due to the parallel 
plate constraint Within the weak-coupling theory of 
BCS p-wave superconductor, we obtain 



2(1 - t) for t r- 
1 - 0{t^) for t 



1 







(5) 



where t = T/T^. For the pressure of 27.3 bar with the 
superfluid transition temperature Tc = 2.46 mK corre- 
sponding to the experimental set-up by Yamashita et al^, 
we found 



1 + 2.295(1 -t). 



(6) 



where (1 - t) « 1. 

When the parallel plates are rotated with the rotation 
axis perpendicular to the plates and rotation speed Q, 
the density of iV = 1 vortices is given by 



(tto^)-^ = 2n/{pK) 



(7) 



where k ~ 271%/ {2m'iHe) ~ 6.59 x lO^^m^/s is the quan- 
tum of circulation, and p is the number of quanta per 
vortex. 2a corresponds roughly the distance between two 
vortices, if we divide the two-dimensional space in circu- 
lar cells with radius a. Therefore a is related to f2 by 
a = y^pK/ (271^1). The free energy for one singular vortex 
in unit cell is given by 



= ttxnC^KIu 



(8) 



which does not involve the texture energy except $ = 
tan~^{y/x). Here ^ ~ 0.01 fj,m is the coherence length. 
On the other hand, the free energy for free pairs of 




FIG. 1: Phase diagram for SV and free HQV phases as func- 
tions of K and a/^. Singular vortex (SV) phase is stable only 
for small a/^ for all temperatures. 



HQVs with their distance R{< 2a) is obtained as follows. 



F, 



fH 



.+ ^/a^ - i?2/4 



2^ 



R . R 

arcsm — — 

4a \ + i?2/4 



(9) 



To get the above free energy, we used the following forms 
for $ and d. 



$ = 
and 

with 



tan 



x + R/2 



tan 



X - R/2 



xcosa -f y sma 



"=2 



tan 



y J 



X- R/2 

y 



(10) 



(11) 



(12) 



As in the case of the singular vortex, this solution has 
dz = 0. If we neglect the last term in the Eq. (9), FfH 
can be easily optimized for R/{2a) = \/2K + 1/{K + 1) 
and 



Ff 



fH 



TT \ a 1 + 2K 



In 



/2a Vl + 2A' 



2Vl + 2K 
l + K 



■ arcsm 



yi + 2K 

V2 + 4/\ + A'2 



(13) 



Now we compare the free energies of singular vortex, Eq 
dHl) and pairs of free HQV, Eq. and find that the 

free HQVs are stable than the singular vortices if a/^ 
satisfies the following equation. 



a ^ 1 + 2/^ 



i - 2(1 + A') V {l + Kf 



{l + 2Kf/'^ \ <^~^> 



(14) 



In the above equation, the last term in Eq. (|13p is ne- 
glected for the analytic solution^^ However, we find that 
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FIG. 2: A configuration of d-vector for /3 = 



Bound HQV 



Free HQV 



this approximation does not make a qualitative difference 
for the phase diagram as a function of (a/^) and K shown 
in Fig 1. The free HQV is stable than the singular vortex 
for all temperatures ( 1 < < 3.295), since we are inter- 
ested in small rotations of << where f2„ ~ 5 x 10^ 
rad/s which corresponds to a ~ ^ ~ 0.01/im. Hereafter, 
we will limit ourselves to HQVs and compare the free 
energies between free HQVs and bound HQVs. 



III. FREE ENERGY OF THE BOUND HQVS 
PAIRS 

If we allow d-vector to bend out of the x — y plane in 
the vicinity of d-digyrations, the system will lower the 
energy. So it is quite possible that two HQVs are bound 
when they are close enough. Let us consider the following 
d configurations. 



d~~ cos a{x cos /3 + y sin /3) + 



sma^ 



(15) 



where P is the cyclic variable, and an example of d con- 
figurations is shown in Fig. 2 for /3 = 0. 

Using the free energy formula for the bound HQVs, 
Eq. (8) in [H, we find 



F, 



bH 



-XnC^< ATn 



a + - i?2 /4 



R_ 

4^ 



R 



■ arcsm 



-t- i?2/4 



2e 

R 




The free energy can be optimized for R/a for given K 
and a/S,H assuming that a/^^f >> 1. We found that 
R/a = 0.3 for a/S^n = 10 for all K values, and R/a = 0.15 
for a/S,H = 20 for aU K values. For a/^n = 30, R/a = 0.1 
for all K values. The ratio, R/£^h{= R/a x a/S^n) ^ 3.0 is 
independent of K and a/S^H- It is remarkable that R/^h 
is independent of temperatures. 

Using the optimal value of R/^h, the free energy of 
the bound HQV can be written as 



^^XnC- \ Kin ( ^ 



3^^ 



In ( 3^^ 



arcsm 

4 a 



1 + 



9?^ 



4 U 



*l-7) 



FIG. 3: Phase diagram for bound HQV phase, and free HQV 
phase as functions of K and a/^H- Note that the free HQV 
is only stable for a relatively small a/^n for all K 



We compare the free energies of bound HQV, Eq. ( \T7} 
and free HQV, Eq. ( [T5)) neglecting the correction terms 
of d textures (the terms with arcsin(..)). We find that 
the bound HQVs is stable over the free HQV, if a/£,H 
satisfies the following equation. 

^ > 1(2)^(1 + i^)i+^(l + 2K)-^''+i'> (18) 
3 

The phase diagram for free HQV and the bound HQV 
phases as a function of a/S,H and K is shown in Fig. 
3. We see that the bound HQV is more stable than the 
free HQV when a > 10£,h for all temperatures. In other 
words, the bound pair HQVs are always more stable than 
free HQVs for relatively small rotation, ft < flc, where 
Clf. rad/s, assuming that ~ 5 /im. From these 
results, we analyze the fl dependence of the satellite in- 
tensity (Ig) as follows. 

In the small rotation speed, the bound pairs are most 
stable. So in the equilibrium condition, the satellite in- 
tensity increases linearly with f2. When Q reaches Oc, 
the free HQV pairs become more stable. Many bound 
pairs are converted into free pairs. This looks a case of 
phase separation. In order to clarify this situation, NMR 
experiment with the field gradient in the parallel plates is 
highly desirable. It is also interesting to note that the 
dependence of the satellite intensity Ig looks surprisingly 
similar to —AnM versus B in typical type II supercon- 
ductors as first formulated by AbrikosovJi 



IV. SPIN WAVE SPECTRUM 

As we mentioned above, both the singular vortices and 
free HQVs do not produce any characteristic NMR sig- 
nals, since their presence does not disturb the configu- 
ration of ^ _L d. On the other hand, the bound pairs 
of HQVs are clearly visible by NMR as discussed in 
[Hfietl. We adopt the formulation for the NMR satel- 
lite frequencies and intensities studied in p^ . However, 
in the present case there will be no potential energy for 
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the longitudinal mode, since the rotation of d around 
the z-axis does not change the dipole energy. In other 
words, the bound pairs of HQVs as well as singular vor- 
tices and free HQVs are invisible by longitudinal NMR. 
All of them give the no longitudinal resonance. In con- 
trast, the transverse NMR is of great interest. 

The transverse resonance frequency, Xg is determined 
as follows. Comparing to [l6| . we consider the bound 
HQVs when H \\ z and < S.d- Let us consider a small 
oscillation of d as follows. 

d + Sd = cos {a + g) [cos(/3 -I- f)x + sin (/3 -I- f)y\ 

+ s\i\{a + g)z. (19) 

The change of free energy due to the small oscillation up 
to the quadratic order in / and g is found as follows. 



5F = Y,\VM,? +CH\5d, 



P = cos2(a)|V/|2 + |V5p 
e^2(i_2sin2(a))5^ (20) 



Therefore, the satellite frequencies of A/ and are ob- 
tained by the following eigenvalue equations. 



- C|fV(cos2 (a)V/) = A/cos2(a)/ 
-e|fV2.g-(l-2sin2(a))5 = A,,^, 



(21) 



As we mentioned above, the longitudinal mode has a triv- 
ial solution of / = constant with A/ = 0. Therefore, 
there will be no longitudinal resonance. As to the trans- 
verse mode, the potential Vg{x,y) = — 1 + 2sin^ (a) is 
shown in Fig. 4, and sin^ (a) is given by 



V„ 




FIG. 4: Potential for the transverse mode is shown in x — y 
plane. It is repulsive between two HQVs, where R is the 
optimal separation between the bound HQVs. 




FIG. 5: The satellite frequency of the transverse mode as a 
function of /R is shown for the bound HQVs. 



sin^ (a) 



■ 2 

sm V 



sinh^ u + sin^ v 



using the change of parameters as 



X = — cosh u cos v, 
2 



y 



R 



sinhiisinw. 



(22) 



(23) 



Note that the potential energy in the present case is 
completely different from the one considered in . This 
implies (a) first of all spin waves are extended, and (b) the 
related intensities are ~ ■na'^n^. At the critical rotation 
of 51 = rJc, ny = (tto^)"^. 

These eigenequations can be solved using a varia- 
tional function. Note that there is a trivial solution of 
50 = sechw with Ag = —1, which does not give a dis- 
tint satellite signal. So we look for a solution which is 
orthogonal to go and have a node in the spin wave func- 
tion. We find the following function of g with a vartional 
parameter ^ satisfies the above criteria. 



(sechw)'^ — p X tanhu x (sechit)^ 



(24) 



where v = 1.13 and p = (^^ _ i)(^b[\,^] - 
\B[\, ^])/(m + 3), and B is the Euler beta function. 
With the above variational functions, we find the opti- 
mal value of Ag which is obtained by minimizing Ag with 



respect to /i. Ag vs. £,h/R is shown in Fig. 5. Note 
that R/^H S is independent of a/^n and K, as long as 
the system is in the regime where the bound HQVs are 
stable. The value of Ag ^ —0.05 for R/£,h = 3 is close 
to the satellite frequency, Ag = —0.06 ~ —0.09, observed 
in [ij. On the other hand, we find that the bound state 
does not exist, when £^h/R> 0.35. 

The intensity of the transvermode is obtained as 



I r 

I Jo 



arccosh{2a / R) 



e'^"gdu\ 



/c 



arccosh{2a/ B.) 




gl'^du 



(25) 



For R/^H = 3, Ig depends on a/£_H- We find that Ig = 
23.4(i?/2)2n„ and Ig = 43.2(i?/2)2n„ for a/£,H = 7, and 
a/S,H = 10, respectively. Since n^ = (tto^)"^ at = fic, 
Ig at the critical rotation is given by 



43.2 / 3 
20 



0.31 for a/^H = 10, 



93 4 / 3 \ 

= —[YaJ = ^-^^ "/^^ = ^- 

The ratio between the satellite intensity and total inten- 
sity, SI /I in the NMR spectra was observed as SI /I = 0.2 
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at T = 0.8Tc, and 0.1 at T = 0.93TcJ. Therefore, the ob- 
served signal is about 60 % and 30 % of the theoretical 
value at T = 0.8Tc and T = 0.93Tc, respectively. In 
our analysis, Ig does not explictily depend on tempera- 
ture, while it is possible that the vortex pinning changes 
the satellite intensity. Therefore, the temperature depen- 
dence in the experiment suggests that the vortex pinning 
becomes weaker as temperature decreses. 

V. DISCUSSION AND SUMMARY 

Since Salomaa-Volovik's proposal for the HQVs in 
1985, HQVs have been searched in the rotating super- 
fluid •^He-A in the confined geometry. However, all of 
earlier experiments have been done in the wider gap D 
(> 25fim). Then the most common texture is the vor- 
tex sheet associated with the composite solitonJ^ In the 
perspective, the vortices seen in Walmsley et aiM ap- 
pears to be the vortex sheet as well. The experimen- 
tal data in [l[ is quite different from the earlier ones in 
the following ways, (a) A new transverse satellite with 

Xg = -0.06 0.09 is observed, (b) Both the satellite 

frequency and the critical value of flc where the maxi- 
mum value of SI /I is observed are almost independent of 
temperature, (c) The critical value of Qi where the first 
appearance of the satellite singal occurs and the maximal 
value of SI / 1 appear to weakly depend on temperature. 

The above observations are in fact consistent with our 
HQV scenario. We have shown that R/^h ~ 3.0 is inde- 



pendent of H and T. Also from the universality of this ra- 
tio, we predict that the satellite frequency and the inten- 
sity are almost independent of H and T, when we ignore 
the vortex pinning effect. On the other hand, any texture 
involving Z-texture should depend strongly on tempera- 
ture and gap size D. Furthermore, another satellite fre- 
quencies associated with Z-bend textures should appear 
near Xg ~ —li^ Absence of these satellites clearly indicate 
that the satellites seen in Yamashita et al^ should be due 
to the bound pairs of HQVs. The parallel experiments 
with different size of D arc highly desirable to distinguish 
the d-texture vs. Z-textures as well, as our study indicates 
that the satellite frequency is independent of D. 
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